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Stirring  and  mixing:  What  are  the  rate-controlling 
processes? 

Chris  Garrett 

Department  of  Physics  and  Astronomy,  University  of  Victoria,  Victoria, 

British  Columbia,  Canada 


Abstract.  The  parameterization  of  unresolved  processes  in  oceanographic 
models  is  complicated  by  the  interplay  of  processes  on  a wide  variety  of  space 
and  time  scales.  The  lack  of  spectral  gaps  further  complicates  the  situation, 
though  perhaps  not  as  seriously  as  might  be  feared.  Understanding  the 
interaction  of  different  processes,  and  determining  which  one  is  critical,  or 
rate-controlling,  may  be  aided  by  thinking  in  terms  of  a triple  decomposition 
into  mean,  eddies  and  turbulence.  Particular  physical  processes  reviewed 
include  the  fate  of  energy  released  in  baroclinic  instability  and  the  ultimate 
thermohaline  frontolytic  mechanism  on  mean  isopycnal  surfaces.  Both  of 
these  issues  may  have  implications  for  diapycnal  mixing  rates.  Another 
important  question  reviewed  briefly  is  that  of  the  efficiency  of  conversion  of 
turbulent  kinetic  energy  to  mean  potential  energy;  we  do  not  know  what 
“external”  parameters  determine  it.  The  kinematic  details  of  overturning  are 
also  discussed  in  terms  of  the  probability  distribution  of  displacements  from 
a stably  stratified  buoyancy  profile. 

1.  Introduction 

The  behavior  of  the  ocean  is  affected  by  a variety  of 
processes  occurring  at  space  and  time  scales  that  are  too 
small  to  be  resolved  explicitly  in  models.  The  influence 
of  these  processes  then  needs  to  be  parameterized  in 
terms  of  variables  that  are  explicitly  included  in  the 
models. 

Following  a “Reynolds  decomposition”  of  variables 
into  a (slowly  changing)  mean  and  fluctuations,  we 
therefore  seek  a representation  of  the  eddy  fluxes  u'C' 
for  a scalar  C and  u'u'  for  momentum,  where  the  prime 
superscript  denotes  the  fluctuation.  For  the  scalar  C it 
is  usual  to  assume  that  the  eddy  flux  is  related  by  a 
tensor  to  the  local  gradient  of  the  mean  C,  as  would 
be  appropriate  in  a mixing  length  theory  in  which  the 
fluctuations  are  related  to  particle  displacements  that 
are  small  compared  with  the  distance  over  which  the 
mean  changes  significantly.  Thus 

= -TijdC/dxj.  (1) 

The  symmetric  part  of  is  diagonalizable  and 
is  likely  to  represent  down-gradient  diffusion  parallel  to 
the  principal  axes  of  the  tensor.  The  antisymmetric 
part  Sij  of  Tij  has  an  associated  “skew  flux”  Fs  given 


by  _ _ 

Fsi  = —SijdC/dxj  = (D  x VC)i  (2) 

where  D = — (S23, 5a  1, Si2)-  This  flux  is  perpendicular 
to  VC  and  may  be  written  as 

Fs  = -(V  x D)C  + V x (DC).  ‘ (3) 

The  second  term  of  this  is  non-divergent  and  so  does 
not  affect  the  evolution  of  C.  The  first  term  is  advective 
with  a velocity  Us  which  may  be  written 

Usi  = dSij/dxj.  (4) 

This  standard  formalism  (e.g.,  Rhines  and  Holland 
1977;  Moffatt  1983;  Middleton  and  Loder  1989)  is  purely 
kinematic.  In  practice,  however,  it  seems  likely  that  Kij 
describes  large  mixing  rates  along  mean  isopycnals  and 
a very  much  smaller  diapycnal  mixing  rate,  whereas  Us 
is  related  to  the  difference  between  Lagrangian  and  Eu- 
lerian  mean  flows.  Values  for  and  Us  appropriate 
for  the  oceans  in  their  present  state  could  conceivably 
be  obtained  from  appropriate  and  sufficient  observa- 
tions. However,  extrapolation  of  limited  measurements 
and  the  requirements  of  models  that  seek  to  be  pre- 
dictive for  different  ocean  states  require  that  we  obtain 
formulae  for  Kij  and  Us  in  terms  of  resolved  variables. 
This,  in  turn,  requires  that  we  understand  the  processes 
responsible  for  the  fluxes. 
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Many  problems  immediately  arise.  One  is  the  va- 
lidity of  a representation  of  the  fluxes  in  terms  of  local 
mean  gradients.  There  are  clearly  examples  in  the  ocean 
in  which  ocean  properties  are  advected  by  coherent  ed- 
dies over  larger  distances  than  those  over  which  mean 
gradients  are  reasonably  constant.  This  is  likely  to  be 
even  more  of  a problem  for  momentum,  which  may  be 
carried  long  distances  by  waves.  Such  situations  call 
for  different  forms  of  parameterization  but  will  not  be 
considered  here. 

Another  fundamental  difficulty  is  associated  with  the 
decomposition  into  mean  and  fluctuations.  This  es- 
sentially assumes  that  there  is  a spectral  gap  of  some 
sort  between  the  mean  and  fluctuations,  allowing  an 
assumption  that  the  mean  changes  little  over  a time  or 
space  scale  large  enough  for  the  determination  of  statis- 
tically accurate  eddy  fluxes.  This  important  issue  will 
be  discussed,  albeit  naively,  in  Section  2. 

A further  basic  issue  concerns  the  extent  to  which 
different  ocean  processes  are  in  parallel  or  in  series. 
If  in  parallel,  which  is  dominant?  If  in  series,  which 
one  controls  the  eddy  flux,  leaving  the  other  just  to 
do  what  it  has  to?  This  will  be  discussed  in  Section 
3,  followed  by  a discussion  in  Section  4 of  the  possi- 
ble insights  to  be  gained  from  a “triple  decomposition” 
into  mean,  eddies  and  turbulence.  Section  5 addresses 
the  inter-relationship  of  isopycnal  stirring  and  diapyc- 
nal  mixing,  and  Section  6 speculates  on  the  factors  that 
may  determine  the  “mixing  efficiency”  with  which  me- 
chanical energy  input,  such  as  that  from  breaking  inter- 
nal waves,  is  converted  to  an  increased  potential  energy 
of  the  mean  state.  Some  clues  on  this  may  come  in  the 
future  from  careful  analysis  of  overturning  motions  in 
the  ocean,  so  some  recent  ideas  and  results  on  this  topic 
are  summarized  in  Section  7. 

2.  The  Spectral  Gap 

At  mid-latitudes,  internal  waves  may  overlap  in  spa- 
tial scale  with  other,  quasi-geostrophic,  motions,  but 
have  different  time  scales.  Thus,  in  determining  the  ef- 
fect of  the  former  on  the  latter,  one  could  presumably 
exploit  the  existence  of  a spectral  gap  in  the  frequency 
domain.  Near  the  equator  the  internal  wave  time  scales 
may  be  longer  and  overlap  with  the  time  scales  of  other 
motions,  so  that  a separation  does  not  appear  to  be 
possible  in  either  time  or  space.  Nonetheless,  a “dy- 
namical gap”  still  exists,  which  is  what  permits  one  to 
talk  separately  about  the  two  classes  of  motion,  and  one 
assumes  that  it  will  still  be  possible  to  parameterize  the 
effects  of  the  waves  on  the  other  motions. 

Motions  other  than  internal  waves  (and  even  smaller 
phenomena  such  as  those  associated  with  double  diffu- 
sive processes)  may,  of  course,  have  a continuous  spec- 


trum, without  any  spectral  gap.  This  would  appear  to 
make  parameterization  of  small  scales  impossible.  The 
success  of  Large  Eddy  Simulations  (e.g.,  Metais  1998) 
has  shown,  however,  that  the  details  of  unresolved  scales 
may  be  unimportant  if  their  only  role  is  to  absorb  vari- 
ance that  is  generated  at  larger  scales  and  cascades  to 
small  scales.  This  could  occur  in  three-dimensional  tur- 
bulence, with  the  only  requirement  being  that  the  start 
of  the  inertial  subrange  be  resolved. 

In  the  ocean,  of  course,  as  in  the  atmosphere,  the  cas- 
cade may  be  partly  to  larger,  rather  than  smaller,  scales 
so  that  the  details  of  the  small-scale  behavior  do  affect 
large  scales  after  a finite  time.  This  is  the  classic  prob- 
lem of  chaotic  behavior  and  lack  of  predictability.  On 
the  other  hand,  experience  in  the  atmosphere  suggests 
that,  while  “weather”  is  unpredictable,  “climate”  may 
be  predictable  (e.g.,  Mote  and  O’Neill  2000).  Maybe 
we  can  hope  for  the  same  in  the  ocean. 

A naive  view  might  thus  be  that  the  lack  of  a spec- 
tral gap  is  less  important  than  might  have  been  feared. 
Either  separation  is  possible  using  a dynamical  gap,  or 
unresolved  scales  do  not  have  a back  effect  on  resolved 
scales,  or  the  back  effect  does  not  affect  the  climatic 
state  of  the  ocean.  The  issue  is,  however,  one  that 
should  be  kept  under  review. 

3.  Are  Processes  in  Parallel  or  in 
Series?  What  is  Rate-Controlling? 

In  three-dimensional  turbulence  at  high  Reynolds 
number,  the  turbulent  scalar  flux  u C'  is  in  parallel 
with  a much  smaller  molecular  flux  — kVC,  with  k the 
molecular  diffusivity.  The  turbulent  flux  does  depend 
on  the  presence  of  the  molecular  diffusivity  in  series;  as 
stressed  by  Nakamura  (1996)  and  Winters  and  D’Asaro 
(1996),  the  total  flux  may  be  thought  of  as  being  purely 
diffusive,  but  across  a highly  convoluted  surface  of  con- 
stant concentration.  The  value  of  the  molecular  dif- 
fusivity does  not  determine  the  value  of  the  turbulent 
flux;  reducing  the  molecular  diffusivity  would  just  lead 
to  an  increase  in  the  streakiness,  or  fine-scale  gradients, 
of  the  scalar. 

Mathematically  this  can  be  summarized  from  the 
equation  for  the  rate  of  change  of  the  concentration 
variance 

anrz  , , 

+V-(uC'2  + u'C"2  - kVC'2]  + 2u'C"  • VC 

= —2k VC'  • VC'.  (5) 

This  is  just  the  Osbom-Cox  formula  ( Osborn  and  Cox 
1972)  if  C is  the  temperature.  If  the  time-dependent 
and  divergence  terms  on  the  left  hand  side  of  (5)  are 
neglected  on  the  grounds  of  stationarity  and  spatial  ho- 
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mogeneity,  there  is  a balance  between  variance  produc- 
tion and  dissipation.  The  turbulent  eddy  flux  may  re- 
main the  same  when  k is  reduced  provided  that  the 
mean  square  gradient  of  the  concentration  fluctuations 
is  increased. 

One  might  say  that  the  turbulent  and  molecular 
fluxes  act  in  parallel,  with  the  former  dominating,  but 
that  the  processes  also  act  in  series,  with  the  turbulent 
stirring  being  the  rate-controlling  process.  The  molecu- 
lar diffusion  just  does  what  it  has  to;  it  is  essential  but 
does  not  control  the  flux. 

One  can  think  of  similar  situations  in  the  environ- 
ment. For  example,  the  meridional  flux  of  potential 
temperature  in  the  atmosphere  is  controlled  by  the  eddy 
flux  associated  with  weather  systems.  The  eddy  flux  re- 
lies on  the  presence  of  radiative  air  mass  modification 
in  series,  but  the  magnitude  of  the  flux  does  not  depend 
on  the  details  of  this. 

On  the  other  hand,  the  heat  flux  across  the  north 
wall  of  the  Gulf  Stream  may  appear  to  be  just  a func- 
tion of  the  rate  of  meander  formation  and  subsequent 
generation  of  warm  core  rings,  which  then  lose  their 
excess  heat  through  air-sea  interaction.  Many  of  the 
rings,  however,  recirculate  and  rejoin  the  Gulf  Stream 
before  they  have  fully  decayed  (e.g.,  Olson  1991).  Pre- 
sumably, therefore,  the  northward  heat  flux  in  this  case 
depends,  at  least  partly,  on  the  actual  strength  of  the 
air-sea  interaction  which  acts  in  series  with  the  initial 
eddy  flux. 


4.  Triple  Decomposition 

When  there  is  some  small-scale  turbulence  acting  in 
series  with  larger  scale  eddies,  some  insight  may  be 
obtained  from  a triple  decomposition  of  the  concen- 
tration field  into  mean,  eddies  and  turbulence  ( Joyce 
1977;  Davis  1994).  It  is  assumed  that  there  is  a spec- 
tral gap  between  the  mean  and  the  eddies,  and  another 
gap  between  the  eddies  and  the  turbulence.  We  write 
C = Cm  + Ce  + Ct,  with  subscripts  to,  e and  t denot- 
ing mean,  eddies  and  turbulence,  respectively.  For  any 
quantity  Q we  write  Q for  the  average  of  Q over  a time, 
or  space  scale,  large  compared  with  the  scale  of  the  tur- 
bulence but  short  compared  with  that  of  the  eddies,  and 
(Q)  for  the  average  of  Q over  a time  or  space  scale  long 
compared  with  that  of  the  eddies  but  short  compared 
with  that  of  the  mean  state. 

Assuming  stationarity,  so  that  variances  do  not  change 
with  time,  and  homogeneity,  so  that  divergence  terms 
may  be  ignored,  it  is  straightforward  to  derive 

(i ^Ct > • VCm  + (ueCe)  • \ (x)  (6) 


Mean 


Eddies 


I 


Turbulence 


Dissipation 


Figure  1.  A simple  schematic  of  the  path  of  scalar  con- 
centration variance  from  mean  to  dissipation  via  eddies  and 
turbulence. 


(u tCt)  ■ NCm  + (u tCt  • VGe)  = - \ ( X ) (7) 

<ueCe)  • VCm  - &ct  • VCe)  = 0 (8) 


where  x = 2k VC"  - VC'  is  the  rate  of  dissipation  of 
scalar  variance,  as  before. 

The  first  equation  here  is  just  the  Osbom-Cox  equa- 
tion if  one  does  no  separate  averaging  over  the  time 
scales,  essentially  lumping  eddy  and  turbulent  fluctu- 
ations together.  The  second  equation,  however,  shows 
that  the  dissipation  may  be  regarded  as  coming  from 
the  production  of  variance  by  the  turbulence  acting  on 
both  the  mean  state  and  the  fluctuations  produced  by 
the  eddies.  The  third  equation  is  just  the  difference 
of  the  previous  two  and  shows  that  variance  generated 
by  the  eddies  acting  on  the  mean  state  must  be  passed 
on  to  dissipation  by  the  turbulence  acting  on  the  eddy 
fluctuations.  The  situation  is  summarized  in  Figure  1, 
showing  the  pathway  from  production  to  dissipation  via 
the  eddies  as  well  as  directly  from  the  turbulence  acting 
on  the  mean  state. 

The  molecular,  turbulent  and  eddy  fluxes  are  all  in 
parallel  in  determining  the  total  flux  at  a fixed  point, 
but  it  is  not  necessary  for  the  eddy  flux  to  dominate 
the  turbulent  flux  in  the  same  way  as  the  turbulent 
flux  dominates  the  molecular  flux  in  the  simple  case. 
Just  as  considered  by  Nakamura  (1996)  and  Winters 
and  D’Asaro  (1996)  for  the  simple  case,  one  may  think 
of  the  flux  in  different  ways.  At  a fixed  point,  or  across 
the  contours  of  the  overall  mean  scalar,  the  flux  is  made 
up  of  the  molecular,  turbulent  and  eddy  fluxes  all  act- 
ing in  parallel.  Alternatively,  the  flux  is  purely  molecu- 
lar across  the  contours  of  the  instaneous  concentration 
field.  A third  way  of  viewing  the  problem,  however, 
is  with  respect  to  contours  of  C obtained  by  averaging 
over  the  turbulence  but  not  the  eddies.  In  this  case  the 
flux  is  carried  by  the  turbulence  and  molecular  diffusion 
in  parallel.  If  the  contours  defined  this  way  axe  not  very 
convoluted  it  implies  that  the  transport  across  a fixed 
surface  is  dominated  by  the  turbulence  rather  than  by 
the  eddies. 

There  are  likely  to  be  variations  on  this  theme,  but 
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it  does  seem  that  a triple  decomposition  may  provide 
insight  in  some  situations,  bearing  in  mind,  however, 
that  there  may  not  be  convenient  spectral  gaps  between 
the  mean,  eddies  and  turbulence. 

5.  The  Relationship  Between  Isopycnal 
Stirring  and  Diapycnal  Mixing? 

We  tend  to  think  of  adiabatic  processes  associated 
with  mesoscale  eddies  in  the  ocean  as  being  indepen- 
dent of  small-scale  diabatic  mixing  processes  caused  by 
things  like  breaking  internal  waves.  It  is  important  to 
consider  whether  this  is  true,  and,  if  not,  what  implica- 
tions it  might  have  for  the  parameterization  of  diapy- 
cnal mixing.  Two  previously  studied  problems  will  be 
briefly  reviewed  here.  The  first  concerns  the  ultimate 
fate  of  the  available  potential  energy  released  into  eddies 
in  baroclinic  instability;  the  second  concerns  the  final 
dissipation  mechanism  for  thermohaline  fronts  that  are 
presumably  generated  by  adiabatic  stirring  on  isopycnal 
surfaces. 

5.1  Does  baroclinic  instability  lead  to  diapycnal 
mixing? 

A popular  parameterization  scheme  for  the  relax- 
ation of  mean  isopycnals  by  baroclinic  instability  in  the 
ocean  was  proposed  by  Gent  and  McWilliams  (1990).  It 
represents  the  horizontal  components  of  the  skew  veloc- 
ity by  ( KVhb/bz)z , where  b is  the  buoyancy  field.  Gent 
et  al.  (1995)  showed  that  the  effect  of  this  is  close  to 
that  of  having  a vertical  eddy  viscosity  Av  = (f2/N2)n 
acting  on  the  mean  flow.  The  appropriateness  of  this 
scheme  over  one  invoking  the  mixing  of  potential  vor- 
ticity  is  discussed  elsewhere  in  this  volume.  The  point 
to  be  made  here  is  that,  as  discussed  by  Tandon  and 
Garrett  (1996),  one  needs  to  consider  the  ultimate  fate 
of  the  mean  available  potential  energy  released  to  the 
eddies.  It  seems  unlikely  that  it  is  dissipated  adiabati- 
cally  by  internal  friction  in  the  ocean.  Some  of  the  eddy 
damping  may  occur  via  air-sea  interaction  or  by  viscous 
damping  in  the  boundary  layer  at  the  sea  floor,  but  one 
also  needs  to  consider  the  possibility  that  the  energy  is 
lost  in  the  fluid  interior  and  that  some  fraction  of  this 
appears,  through  diapycnal  mixing,  as  an  increase  of 
the  basic  potential  energy. 

As  discussed  by  Tandon  and  Garrett  (1996),  the  rate 
of  release  from  the  mean  flow  is  ( f2/N2)K,\uhz\ 2 for  an 
equivalent  eddy  viscosity,  as  above,  acting  on  the  verti- 
cal shear  of  the  mean  horizontal  flow  u?,.  Using  the 
thermal  wind  equation  |u/tZ|  = 1 1 this  energy 

loss  rate  becomes  s2kN2  where  s = | V^bj/iV2  is  the 
mean  isopycnal  slope.  If  the  rate  of  creation  of  mean 
potential  energy  is  at  a rate  Rf  times  this,  and  is  ex- 


pressed as  KVN2,  then  the  vertical  eddy  diffusivity  is 
Kv  = RfS2n  and  could  be  as  large  as  10~3m2s-1  in 
places  like  the  Southern  Ocean.  Tandon  and  Garrett 
(1996)  argue  that  this  could  have  a significant  effect  on 
tracers,  though  it  is  unlikely  to  significantly  augment 
the  spin-down  of  the  mean  flow  already  associated  di- 
rectly with  the  Gent  and  McWilliams  mechanism. 

The  problem  thus  seems  to  warrant  continued  inves- 
tigation. Strong  internal  wave  activity  and  inferred  high 
dissipation  have  been  reported  for  the  Southern  Ocean 
by  Polzin  and  Firing  (1997).  They  suggest  that  the 
waves  originate  as  lee  waves  generated  by  strong  flows 
over  a rough  sea  floor,  and  that  these  strong  flows  in- 
clude currents  associated  with  the  mesoscale  eddies.  In 
that  case  the  internal  waves  would  indeed  be  associated 
with  the  decay  of  the  baroclinic  eddies,  though  the  mix- 
ing in  the  water  column  would  not  be  associated  with 
the  local  current  shear.  It  is  worth  considering  whether 
there  could  additionally  be  some  direct  local  connec- 
tion between  the  mesoscale  eddy  current  shear  and  the 
excitation  of  internal  waves. 

5.2  Are  thermohaline  fronts  dissipated 
passively  or  actively? 

Compensating  lateral  gradients  of  potential  temper- 
ature and  salinity  can  exist  on  isopycnal  surfaces  in  the 
ocean.  Purely  adiabatic  stirring  of  these  gradients  will 
presumably  lead  to  density-compensated  thermohaline 
fronts.  Some  small-scale  mixing  mechanisms  must  oc- 
cur to  prevent  the  temperature  and  salinity  gradients 
from  increasing  without  limit. 

The  gradients  could,  of  course,  become  sharp  enough 
that  molecular  diffusion  alone  removes  them.  Alterna- 
tively, there  might  be  some  small-scale  lateral  mixing 
mechanism,  perhaps  involving  vortical  modes,  which 
would  hasten  the  transfer  to  molecular  scales.  Both 
of  these  possibilities  could  be  described  as  “passive”  in 
the  sense  that  the  small-scale  mixing  mechanism  was 
in  place  already,  and  just  copes  with  the  extra  variance 
produced  by  the  isopycnal  stirring. 

Another,  more  likely,  passive  mechanism  was  ex- 
plored by  Haynes  and  Anglade  (1997)  and  is  reviewed 
by  Haynes  in  this  volume.  The  fundamental  idea  is 
that  the  vertical  shear  of  the  stirring  process  leads  to 
the  development  of  vertical  scalar  gradients  at  the  same 
time  as  lateral  gradients.  Sharp  gradients  develop  in  the 
vertical  as  well  as  laterally,  with  diapycnal  mixing  then 
acting  to  smooth  them.  Haynes  and  Anglade  (1997)  ar- 
gue that  a lateral  frontal  width  of  a few  kilometres  is  a 
plausible  outcome  for  reasonable  values  of  the  vertical 
shear. 

An  alternative  mechanism  for  the  thermohaline  ffon- 
tolysis  was  proposed  by  Garrett  (1982)  and  reviewed 
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by  Garrett  (1989).  He  suggested  that  as  the  front  de- 
velops, it  will  be  unstable  to  thermohaline  intrusions 
which  would  limit  the  further  narrowing  of  the  front, 
effectively  dissipating  the  lateral  variance  production 
by  larger-scale  stirring.  This  might  be  termed  an  “ac- 
tive” process  in  that  it  introduces  extra  diapycnal  mix- 
ing that  would  not  have  occurred  without  the  lateral 
stirring.  The  associated  diapycnal  eddy  diffusivity  is 
inevitably  negative  for  density  (as  is  the  case  for  den- 
sity for  any  process  driven  double-diffusively),  but  is 
likely  to  be  positive  for  salinity  and  can  take  either  sign 
for  temperature,  depending  on  the  circumstances.  The 
process  deserves  further  study  if  the  effective  diapycnal 
diffusivities  that  one  would  add  to  a model  for  large- 
scale  oceanic  behavior  reach  significant  values.  If  the 
diffusivities  are  very  small  one  would  argue  that  the 
process  is  interesting  scientifically  but  just  does  what  it 
has  to  without  having  any  large-scale  impact. 

The  formula  suggested  for  the  effective  diapycnal  dif- 
fusivity of  salt  is 

Ks  = 10~3  D2N(g0\NhS\/N2)3  (9) 

where  D is  the  diameter  of  the  mesoscale  eddies  doing 
the  stirring,  /3  the  density  coefficient  for  salinity  and 
|V/jS|  is  the  lateral  gradient  of  the  large-scale  salinity 
field  (over  a scale  larger  than  that  of  the  eddies  causing 
the  stirring).  The  value  of  Ks  is  clearly  very  sensitive 
to  the  lateral  salinity  gradient  and  only  likely  to  be 
significant  in  a few  places  such  as  the  Mediterranean 
salt  tongue  where  it  is  large  ( Garrett , 1989). 

The  model  also  predicts  the  thickness  of  the  intru- 
sions to  be  about  ^(gp\NhS\/N2)D.  This  could  be 
compared  with  data  in  any  systematic  examination  of 
CTD  profiles.  Of  course  the  intrusions  would  only  be 
found  in  the  narrow  frontal  regions.  The  width  of  these 
was  predicted  to  be 

W ~ 0.08(5£|VfcS|/AT2).DjVfr 1 (10) 

where  fl  is  the  large  scale  strain  rate  (say  10-6s-1  or 
less  typically). 

Presumably  this  active  frontolytic  mechanism  only 
occurs  if  the  width  from  (10)  is  larger  than  the  width 
in  the  passive  mechanism  proposed  by  Haynes  and 
Anglade  (1997).  The  particular  values  for  W tabu- 
lated by  Garrett  (1989)  show  that  this  may  be  the  case, 
though  clearly  there  is  considerable  uncertainty  in  both 
models. 

One  location  for  which  information  is  available  is  the 
site  of  the  North  Atlantic  Tracer  Release  Experiment 
( Ledwell  et  al.  1993).  In  this  experiment  an  artificial 
tracer  injected  at  a depth  of  about  300  m was  teased  out 
into  streaks  by  a strain  rate  estimated  as  3 x 10-7s_1 
and  reached  a width  of  about  3 km.  This  width  was 


found  by  Haynes  and  Anglade  (1997)  to  be  plausible 
for  the  mechanism  they  described,  though  with  consid- 
erable uncertainty  associated  with  the  lack  of  informa- 
tion on  the  vertical  shear  of  eddy  currents. 

To  apply  the  model  considered  here,  we  need  to 
assume  that  at  the  same  time  as  the  tracer  became 
streaky,  the  associated  lateral  convergence  would  lead 
to  thermohaline  frontogenesis.  The  large-scale  isopy- 
cnal  salinity  gradient  in  this  region  is  approximately 
10"7m_1,  and  AT2  ~ 1.8  x 10"5  so  that  g(3\NhS\/N2  ~ 
4 x 10-5.  If  we  take  D ~ 100  km,  then  (10)  gives  a 
frontal  width  of  about  5 km.  This  is  also  close  to  the 
observed  width.  A point  against  this  interpretation  is 
that  no  intrusive  features  were  reported,  though  Gar- 
rett’s (1982)  model  would  predict  them  to  be  only  2 m in 
vertical  extent  with  very  small  temperature  and  salinity 
signatures  and  so  perhaps  not  readily  observable. 

It  also  seems  quite  likely  that  in  this  particular  sit- 
uation neither  of  the  models  discussed  here  constituted 
the  ultimate  lateral  mixing  process,  but  that  the  mix- 
ing was  actually  accomplished  by  small-scale  vortical 
modes.  The  discussion  is  perhaps  somewhat  academic 
anyway;  even  if  Garrett’s  (1982)  mechanism  were  the 
ultimate  frontolytic  process  in  the  NATRE  region,  the 
regional  average  diapycnal  mixing  rate  for  salt,  using 
(9),  would  be  less  than  10~8m2s-1  and  so  utterly  neg- 
ligible. (The  smallness  of  this  value  is  partly  associ- 
ated with  the  small  area  occupied  by  the  thermohaline 
fronts.) 

Perhaps  the  realistic  conclusion  at  this  stage  is  that 
we  do  not  know  in  general  what  the  frontolytic  mecha- 
nism is,  that  it  may  vary  from  place  to  place,  and  that 
in  a few  situations  it  may  be  important  to  know  what 
it  is  if  it  introduces  further  diapycnal  mixing. 

6.  Mixing  Efficiency 

An  important  issue  in  considering  the  relationship 
of  stirring  to  mixing  in  the  diapycnal  sense  is  that  of 
the  mixing  efficiency:  how  much  of  the  energy  being 
put  into  stirring  motions  is  converted  to  the  potential 
energy  of  the  basic  state?  The  topic  is  discussed  by 
Barry  elsewhere  in  these  proceedings  and  also  by  St. 
Laurent  and  Garrett  (2001). 

The  standard  approach  is  to  write 

Kv  = Te/N2,  T = Rf/(  1 - Rf),  (11) 

where  Rf  is  the  flux  Richardson  number,  representing 
the  fraction  of  energy  put  into  the  stirring  motions  that 
produces  a diapycnal  buoyancy  flux  rather  than  being 
dissipated. 

It  is  customary  to  choose  T to  be  about  0.2,  based  on 
comparisons  of  different  methods  of  measuring  the  di- 
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apycnal  diffusivity  (e.g.,  Oakey  1982),  but  there  seems 
to  be  no  reason  why  it  should  be  a universal  constant. 
Ivey  and  Imberger  (1991)  argue  that  F varies  consider- 
ably in  a two-dimensional  space  described  by  two  di- 
mensionless parameters.  One  is  the  “overturn  Froude 
number”  Ftt  = (Lr/Lq)2^3,  where  Lr  = (e/TV3)1/2 
is  the  Ozmidov  scale  and  Lc  is  “the  scale  of  the  most 
energetic  overturns”.  The  other  parameter  is  the  “over- 
turn Reynolds  number”  Rer  = ( Lc/LrY /3,  where 
Lk  = (i^/e)1/4  is  the  Kolmogorov  scale. 

We  note,  however,  that  Fry  is  an  “internal”  dimen- 
sionless parameter,  with  its  value  determined  by  the 
turbulent  mixing  itself  rather  than  by  the  physical  pro- 
cesses driving  the  mixing.  It  thus  does  not  seem  to  be 
a satisfactory  descriptor.  We  might  also  expect  that  F 
should  be  independent  of  Rer  for  large  values  of  the  lat- 
ter or,  equivalently,  of  the  turbulent  Reynolds  number 
e/(vN2)  — Fr^Rer  when  this  is  large. 

Thus  one  still  seeks  “external”  dimensionless  param- 
eters on  which  T might  depend.  One  wonders  whether  F 
might  be  less  if  the  shear  is  persistent  (compared  with 
IV-1)  at  a particular  location;  then  for  much  of  the 
time  the  turbulence  is  mixing  water  that  is  already  well 
mixed,  hence  with  less  buoyancy  flux  resulting.  This 
persistence  might  be  associated  with  the  frequency  con- 
tent of  the  internal  wave  spectrum  (with  an  enhance- 
ment of  low  frequency  energy,  such  as  that  at  inertial 
and  tidal  frequencies)  increasing  the  persistence  and  re- 
ducing T.  Alternatively,  persistence  could  be  a conse- 
quence of  greater  exceedance  of  some  stability  criterion: 
if  the  Richardson  number  goes  well  below  1 /4  at  its  min- 
imum, then  its  duration  below  1/4  is  increased  over  a 
case  in  which  it  just  barely  goes  below  1/4. 

Perhaps  an  apparent  dependence  of  T on  e/(i/N2) 
might  just  be  because  this  parameter  is  acting  as  a 
proxy  for  some  other  indicator  of  mixing  strength  that 
does  not  involve  the  viscosity  u. 

A test  in  the  field  could  of  course  come  from  simul- 
taneous data  on  e and  the  temperature  microstructure 
dissipation  rate  x if  both  were  available,  as  in  Oakey 
(1982).  If  only  velocity  microstructure  data  and  coarser 
density  data  were  available  one  could  proceed  via  a com- 
parison of  Thorpe  and  Ozmidov  scales.  The  former,  Lt, 
is  the  r.m.s.  displacement  in  reordering  a density  profile 
to  be  stable.  The  latter,  Lo,  is  defined  as  (e/TV3)1/2.  If 
we  write  Lo  = uLt,  then  Kv  = Te/TV2  implies  that 

Kv  = Va2NL2T.  (12) 

With  T ~ 0.2  and  a — 0.79  ( Dillon  1982)  or  0.66  ( Craw- 
ford 1986)  this  gives  Kv  ~ O.lTVLf-.  This  latter  result 
seems  likely  to  be  a consequence  of  the  dynamics  of 
overturning  and  more  likely  to  be  a general  result  than 
a fixed  proportionality  between  Lq  and  Lt-  Hence  if  a 


is  found  to  differ  from  0.8,  one  could  use  (12)  to  imply 
T ~ 0.1a-2.  Interestingly,  Ferron  et  al.  (1998)  found 
a ~ 0.95  in  the  energetic  abyssal  mixing  in  the  Ro- 
manche  Fracture  Zone.  As  they  point  out,  the  error 
bars  on  all  these  values  of  a are  such  that  they  could  be 
consistent  with  a universal  value,  but  we  note  that  the 
higher  value  they  find  might  also  imply  a lower  mixing 
efficiency. 

7.  The  Probability  Distribution  of 
Thorpe  Displacements 

The  existence  of  statically  unstable  portions  of  a wa- 
ter column  is  clear  evidence  for  vertical  stirring  and 
mixing.  The  nature  of  this  stirring  is  usually  charac- 
terized by  the  single  parameter  of  the  Thorpe  scale,  as 
discussed  above.  Recently,  however,  Stansfield  et  al. 
(2001)  have  suggested  that  more  information  might  be 
available  from  the  actual  probability  distribution  Pi  ( L ) 
of  the  Thorpe  displacement  L (the  particle  displace- 
ment in  the  reordering  process).  The  Thorpe  scale  Lt 
is  simply  (/0°°  L2P\(L)dL )1/f2  . Now 

Pi(F)  = f HP3(L/H)P2(H)dH/  [ H2P2(H)dH 
Jo  Jo 

(13) 

where  P2  (H)  is  the  probability  distribution  of  overturns 
of  thickness  H (an  overturn  being  a closed  set  of  dis- 
placements) and  P3(L/H)  is  the  normalised  probability 
distribution  of  displacement  L. 

Stansfield  et  al.  (2001)  do  not  offer  any  model  for 
P2(H),  but  suggest  a very  simple  kinematic  model  for 
P3(L/H),  based  on  the  assumption  that  each  particle 
has  an  equal  probability  of  going  to  any  other  location 
in  the  overturn.  This  is  equivalent  to  the  assumption 
that  any  rearrangement  of  the  elements  that  make  up 
an  overturn  is  equally  likely.  If  the  overturn  is  made  up 
of  n points,  then  we  may  take  L/H  = m for  0 < m < n 
and  it  is  easy  to  show  that  P(0)  = 1/n  and  otherwise 
P(m)  = 2(n  — m)/n2.  Then  Lj,  = H2/6  for  large  n. 

There  are  subtleties  associated  with  this  model  as 
it  is  necessary  to  exclude  the  rearrangements  that  do 
not  constitute  a single  overturn,  but  these  are  not  crit- 
ical if  n is  larger  than  about  20  or  so.  Figure  2 shows 
P3{L/H)  for  winter  and  summer  data  from  Juan  de 
Fuca  Strait.  While  the  error  bars  are  not  shown  here, 
there  are  clearly  more  small  displacements,  and  fewer 
medium  and  large  displacements,  in  the  data  than  in 
the  model.  There  is  little  reason  to  regard  the  model 
as  anything  other  than  a convenient  reference,  but  it 
seems  possible  that  examination  of  P3(L/H)  and  also 
P2  ( H ) in  different  settings  might  reveal  different  mixing 
regimes  and  possibly  provide  an  indicator  of  different 
mixing  efficiencies. 
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Figure  2.  The  probability  distribution  of  displacements, 
normalised  by  the  overturn  height,  from  winter  1998  (thick 
line  ) and  summer  1998  (thin  line)  data  from  Juan  de  Fuca 
Strait.  The  dashed  line  is  the  prediction  of  a simple  kine- 
matic model,  assuming  50  points  in  a profile  (though  this 
number  only  affects  the  value  of  L/H  at  which  the  model 
prediction  rises  to  from  1 to  2.)  From  Stansfield  et  al. 
(2001). 

Conclusions 

This  short  review  has  attempted  to  raise  some  of 
the  issues  of  stirring  and  mixing  in  the  ocean,  bearing 
in  mind  that  a primary  goal  is  to  derive  parameteri- 
zations  of  small-scale  processes  that  have  a significant 
impact  on  important  aspects  of  the  ocean’s  behavior. 
An  overriding  issue  is  the  need  for  formulae,  not  just 
numbers,  for  any  parameterization  in  a model  that  as- 
pires to  predictive  capability.  If,  in  a model  of  a partic- 
ular resolution,  several  sub-grid-scale  processes  appear 
to  be  acting  together,  it  is  clearly  important  to  decide 
which  processes  are  rate-controlling  so  that  theoretical 
and  observational  programs  can  be  focused  on  them. 
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